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Abstract 

Properties of the So superfluid phase are studied for symmetric nu- 
clear matter at finite temperature. It is described within a covariant 
hadronic field model, of the a — oj type, with addition of density depen- 
dent correlations simulating effects due to finite extension of nucleons. 
The model is solved in a selfconsistent Hartree-Bogoliubov approach, as- 
suming instantaneous interactions in the superfluid phase. A comparison 
with the results obtained from several hadronic field models is made. Main 
characteristics of our description of the superfluid gap are in qualitative 
agreement with some studies using microscopic potentials, although fur- 
ther refinements could improve its performance. 

1 Introduction 

Superfluid states in the nuclear environment have been extensively studied as 
they have a significative role in several physical processes, such as the struc- 
ture of nuclei out the stability valley and the cooling dynamics of proto-neutron 
stars. 

A variety of models and approximations have been used for this purpose [TJ [21 
[31 [H [H El El El HI HOI HU H21 H31 HH HH HB], mainly non-relativistic potentials 
or effective forces, such as the density-dependent Skyrme or Gogny ones. It is 
not unusual a mixed treatment that combines both schemes in order to simplify 
involved calculations [TJ [2] . 

Approximately two decades ago a covariant model of the field theory of hadrons, 
generally known as Quantum Hadro-Dynamics (QHD) [T71[T5], was used for the 
first time to study nuclear matter superfluidity [5] . There are several reasons to 
use this theoretical framework to deal with nuclear superfluidity, in first place 
there are practical reasons, some self-consistent calculations are more simply 
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stated and easily solved. This property gave rise to a version of QHD, known 
as Density Dependent Hadron Field Theory [19] which casts Dirac-Brueckner 
outputs in the QHD language. In second place one must consider field theory 
as a more adequate tool to make contact with the fundamental theory of strong 
interactions and the fact that covariant formulations are desirable for astrophys- 
ical applications, among others conceptual reasons. Furthermore, the formalism 
has the versatility to include vacuum effects and finite renormalizations in a co- 
herent way [T3J [TS] . 

Since QHD models are formulated as a many-body theory, one of its basic 
premises is the reproduction of the nuclear matter saturation properties. This 
can be fulfilled with a few adjustable parameters and even at the lowest order 
of approximation. Once the free-parameters have been fixed, the QHD model 
has a noticeable prediction power. 

Although a wide spreading of numerical results for the superfluid phase in infi- 
nite nuclear matter can be found in the literature, there is a qualitative agree- 
ment that the superfluid gap should not exceed 3 MeV, and it should vanish 
for densities around the saturation density [13] . Unfortunately, the findings of 
[9 for superfluid nuclear matter in the 1 5'o phase do not agree with these ex- 
pectations. A maximum A max (pp) rs 10 MeV was found there, and only an 
unphysical reduction of 15% in the omega-meson mass yield results comparable 
with currently accepted values. 

This situation has not been changed substantially after the evaluation of dif- 
ferent corrections into the original scheme. However, the good properties of 
this treatment have motivated mixed descriptions combining QHD models and 
conventional potentials [TJ [SJ [3] . 

More recently it was claimed that a coherent inclusion of meson proper self- 
energy [10 , or adjustable quenching factor [IT], could bring numerical calcula- 
tions to the likely values . Further studies about the effect of corrections of the 
meson propagators and the influence of low density instabilities on the super- 
fluid gap can be found for instance in [T2] , 

It was stressed in reference [5] that the high momenta behavior of the repul- 
sive potential has a crucial role in the exceedingly large values obtained for 
^■maxijpF)- It must be noticed that the relative strength of repulsive and at- 
tractive contributions have been calibrated for momenta below the Fermi surface 
in order to produce the saturation mechanism. Therefore it would be desirable 
a pairing potential which preserves the relative strength of its components in 
the Fermi sphere, but having a repulsive component decreasing as faster as the 
attractive one in the high momenta domain. 

Taking these facts into account we try in the present work, to obtain an effec- 
tive and concise model able to deal with the superfluid phase of nuclear matter. 
For this purpose we use and compare several models of relativistic nuclear fields 
interacting through scalar and vector mesons. As a first approach, we reduce 
to its minimal expression the complexity of the nuclear interaction, but further 
refinements can be considered. We introduce a characteristic length scale in 
the effective interaction, which could eventually be traced back to the confining 
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mechanism of the fundamental theory of strong interactions. A similar approach 
was applied in the past to describe heavy ion collisions [20] [21] [22] [23] [24] as 
well as nuclear matter properties [351 023 123 12H1 130] • It has been known 
as finite volume correction since it takes into account the spread of the nucleon 
localization. The relevance of this effect upon the evaluation of some bulk prop- 
erties of the nuclear matter has been stressed long time ago [3T] . 

In the next section we present the theoretical deduction of the general gap 
equation at finite temperature in a context of Landau-Fermi liquid. The last 
part of this section is devoted to the theoretical deduction of the self-consistent 
expression defining the superfluid gap for the 1 Sq phase within the QHD frame- 
work. In section 3 we show the numerical results and discuss them. Finally, 
conclusions are drawn in section 4. 



2 The Formalism 

2.1 Superfluid states in a Landau- Fermi liquid 

The different superfluid phases in a fermionic system can be described in a gen- 
eral and compact way within the formalism of the Fermi liquid, using for instance 
the formalism of reference [32]. It is assumed that the low- lying excitations of 
the system are represented in terms of quasi-particles and, circumstantially, col- 
lective modes. 

We use the notation f a for the equilibrium distribution function of a quasi- 
particle state, where the label a comprise spin, isospin, and momentum quan- 
tum numbers. The fermionic contribution to every conserved quantity, such as 
particle number and energy, can be expressed in terms of a summation over f a . 
In the following we will be interested in nucleon pairs coupled to singlet spin 
and triplet isospin, so that an anomalous distribution function gp and a energy 
gap are introduced. According to [35], we make the decomposition 



gp = g k (p) r( fe M 2 > *™, (1) 



(2) 



, s s' 
i b 



Ap = A k ( P ) r«r( 2 ) *™ t (2) 



(2) 



, s s 
i b 



where momentum, isospin (a, 6), and spin (s, s') dependencies has been clearly 
distinguished. Here Tfc, <7k, k — 1,2,3 stands for the Pauli matrices for isospin 
and spin, respectively. 

On the other hand, the normal phase is filled with quasi-particle states in a 
isospin duplet described similarly by fp — 6 ss >f a b(p), with 

fabip) = fo{p)S a b + fa(p)r^. 
Accordingly, the quasi-particle spectra is assumed in matrix form e(j>) = £q(p) + 

£3(P) 7"3- 
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In the Landau theory of Fermi liquids, the energy of the system E is consid- 
ered as a functional of the distribution functions f a , gp. First variations respect 
to them give the quasi-particle and gap spectra matrix components 



SE 



6 fjk{p)' 



SE 



(3) 



Since / and g are itself functions of e, A, the equations above are self- 
consistent relations. 

Within the block diagonalization procedure of [35] the distribution functions 
are written 



/ = Kn + X(l-n*)X^K, 
g = (KnXy + K(l - n)X, 

He^- XA, ) 1_I 

-l 



K = (l + XX 



(4) 
(5) 

(6) 
(7) 



here t indicates matrix transposition, j3 — 1/kT, £ = £ + /i, and the diagonal ma- 
trix fi =diag(/ii, H2) collects the proton (1) and nucleon (2) chemical potentials. 
The unknown matrix X = XjTjT^o^ satisfies the condition 

£X + X£ T + A - XA^X = 

We have solved this system of equations for symmetric nuclear matter cou- 
pled to T z — 0. We obtained 



Xi = A 

03 (P) = 



; 9i z 
A 3 



= 0, for i= 1,2, X 3 
tanh(/3S A /2), f = 



e - mo ± E A 



h 



2E A 

0, E A = v /A2 + (e - / i ) 



"3 

£o — Mo 
E A 



t&nh(/3E A /2) 



in the first line /io stands for either the proton or the neutron chemical potential. 
The particle density of protons (k=l) or neutrons (k=2) is given by 



rik 



d 3 p 
(2^p 



hip) 



r (3) 
r kk 



MP) 



these equations are used to relate the chemical potentials to the conserved 
isospin and baryonic number density. 

In the next subsection we show the model which provides the quasi-particle 
interaction and spectra. 



2.2 Hadronic models 

Models of the nuclear interaction, formulated in the covariant field theory, have 
been widely used in the study of the dynamics and structure of infinite matter 
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as well as finite nuclei. Since the pioneering work of references [TTJ [TS], the 
simple (7 — lj model has grown in different directions and it was completed in 
order to cover a multitude of manifestations of the nuclear force. 
In particular the subject of the nuclear superfluidity was first treated within 
this context in reference [S] , by using the original a — oj model plus a pseudo- 
scalar pion interaction. The scheme of approximation used there consisted in 
a mean field treatment of the meson and nucleon fields, a Gorkov factorization 
of the pairing interaction and a instantaneous assumption which allows a time- 
independent resolution of the gap equation. Subsequently, this procedure was 
extended to consider the effect of vacuum, and the variation of the in-medium 
meson properties [HI US] . 

Most of these studies agree in a excessively large value for the gap in infinite 
nuclear matter. The realization of a correlated state of two nucleons is a con- 
sequence of the equilibrium between a repulsive and an attractive component 
of the pairing potential, originating in the exchange of virtual w-mesons and 
er-mesons respectively. This mechanism is also found in the binding energy of 
nuclear matter. 

The slow decrease of the repulsive potential as a function of the transferred mo- 
menta, has been pointed out as the main cause of the mismatch. This situation 
can not be modified without a substantial redefinition of the couplings, which 
should lead to a destruction of the saturation mechanism. 

In this work we adopt the simplest version of nucleons interacting through 
a and uj mesons, whose lagrangian density is 

C = * a (i - M a + g s a - g w <Jj) + i (d^ad^a -m 2 s a 2 )- -F^F^ + ^m 2 w iu 2 

where summation over the repeated isospin index a is assumed, F^ v = d^w u — 
dyUi^ and g s ,g w are adimensional coupling constants. The equations of motion 
of the classical fields are 

(i - M a + g s <7 - g w 0) * a = 0, (8) 
(□ + m 2 ) a - g s y a ^ a = (9) 
9 l/ F; + m> M - 5u) * a 7M * a = (10) 

Denoting by D(x,y), D^ v {x,y) the propagators of scalar and vector-meson 
fields, the eqs. © and ([TUf can be formally solved as 

a(x) = g s J d 4 yD(x,y)^ a (y)^ a (y), (11) 

w M (a:) = 9w J a 4 yD^(x,y)y a (y)Y* a (y) (12) 

The Hamiltonian density % is given by the canonical procedure, and the 
energy density of the system £ is evaluated by taking its expectation value 
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£ =< ^ a t lo d ^ a + d ad a-F ^d^ + ^a(a + ml)a-^(d^+mlujn > 

(13) 

Under the hypothesis of static meson fields the second and third terms van- 
ish. 

Inserting the meson equations of motion ©, (|10p together with the Eq. (fl 
into (fT"3"|) . we obtain 

£{x) = < ^{ x ) llo d^ a {x)>+ g l I a i yD(x,y)<*Z(x)* a a (x)* b f) (y)* b f) (y)> 



9& 



d A y D^(xM all ^,) < ^ a a (x)K'(^(y)^Uy) > ( 14 ) 



At this point we introduce an expansion of the nucleon fields in terms of 
quasi-particle creation and annihilation operators, similar to that of a free field 

/M* 
s " p Ejf) i ba M<(p) e ~ iPaX + 4 a (p)<(py PaX ] , (15) 

/M* 
d3p E^pj H a (P)<(Py PaX + d a s ( P mP)e- tPaX ] , (16) 

<{?) = A + M a =u s {Q) 

^2M*(M*+E a (p)) 
P 4- M* 

y/2MZ(M£ + E a (p)) 

but creation and annihilation operators are referred to the lowest energy state of 
correlated nucleons. We have used the quasi-particle properties M* = M — g s s, 
E a (p) = y/M* 2 + p 2 , P a = {E a (p) + g w w,p). Due to the isospin invariance of 
the interaction proton and neutron properties are actually independent. Fur- 
thermore, as we are interested in isospin symmetric matter proton and neutron 
are indistinguishable and the isospin index a becomes superfluous. From now 
on we will omit it, and a degeneracy factor 2 will be included when necessary. 
The quantities s, w stand for the mean values of the a and uj meson fields in 
homogeneous, isotropic matter. As usual, they can be deduced from Eqs. (|9]), 
(jlOp by neglecting derivatives and taking expectation values of the fermionic 
bilinears s — 2 g s < > /m 2 , w — 2 g w < *7 * > /m 2 w . 

As a part of the approximation we neglect in Eq. (|14[) the contribution of 
particle- antiparticle or antiparticlc-antiparticle terms. We define the equilib- 
rium distribution functions for the normal and superfiuid phases 

M* 

- 7mm <bUk)bM> ' m 

M* 

*,(p,*) = 7 ——< b Ak)Up)>, (18) 

M* 

= 7mm <bUmp}> - (19) 



G 



It must been taken into account that f(k,p) oc (27r) 3 (5 3 (p — fc), whereas pairs of 
nucleons are assumed to couple to zero momentum, so that g{k,p) oc (27r) 3 S 3 (p+ 
k). 

Within the mean field approach, corrections to the meson propagation in 
the nuclear environment are dismissed, although a random phase approximation 
could be considered, as in [10]. In the first case, we obtain 



D{z) = / d' l qD{q)e~ lqz (20) 

D(q) = -(qxq^-ml + ie)- 1 , (21) 

D^(q) = -{g^u- q^qv/m 2 w )/{q x q X -m 2 w + ie), (22) 

The term proportional to q^ in (I22p . produce zero contribution in integrals 
combining D M „ with nucleon fields, because of the conservation of the baryonic 
current. 

We are interested in static homogeneous matter, therefore z = x—x' in Eq. ([20)) . 
If we neglect time retardation in the meson propagation, i. e. = xo — x' , then 
terms containing qo are absent in Eqs. (f2"Tj) , (|22p as well as in the exponential 
of the Fourier transform shown lines above. 
All these elements together in Eq. (TUT) produce 

£ = £mfa + £f + £s 



'MFA 



= 4 / a 3 P s(p)f(p) + 2 (Mi J a 3 p-^f( P ) 



2 



,2 I _r3„ f(p)f(l) M* 2 + S(p)S(g) p q 



2 ^ / 



2E(p)E(q) [ e (p) _ £ (g)] 2 - | p - q |2 _ m 2 + fe 



2 /,3 /(?)/(?) 2M* 2 -£(p) J E(< Z ) + p.q 



-2 5 2 / d d g , N 5 — — — £ — 1 (24) 

' ' * H E{p)E{q) [ e (p) _ e ( g )] 2 _ | p_ q |2 _ m 2 + ^ 

■j / r a ^3 gHp)g{q)^ s (p,q) , , 

£s = -g s d a pd a q—- —2 — ■ — -vip,q) 



e(p) - - | p - q | 2 -m 2 

+g 2 Ja 3 pa 3 q MUti^M v{pq) (25) 

J [e(p) - e(g)] - | p - q | 2 -m 2 + ze 

Here we have separated mean field (MFA), Fock (F) and superfluid (S) 
contributions. Certain integrals appearing in (fT4|) vanish because of the isotropy 
of infinite matter. The following notation is used 

X 2 (p,q)~2X(p, g)p-q+pV 



K(p,q) 



E(p)E(q)X(p,q) 

3X(p,p)q 2 + 3X(q, q)p 2 - 2X(p, q)p ■ q + p 2 q 2 



KM = AsiP ' q) ' E { p)E(q)X { p,q) 
r](p, q) = exp (-i2x [e(q) - e(p)}) , 
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with X(p, q) = (M* + E{p)){M* + E(qj). 



The factor p • q = p qv entering in these expressions can be used for integra- 
tion respect to v. For this purpose it is useful the distribution identity 

Denominators in the integrands of Eqs. (|25|) . ([24| come from the meson propa- 
gators, in particular the combinations e{q) — e(p) come from its go dependence. 
In the next step we apply the instantaneous approximation [9 , which result in 
the elimination of all these combinations. 
Within this approach Eq. ([25]) can be re-written 



e 2 [ j3 j3 9(q)gHp) 

£S = - gs J dpdq E(pjW) 



„2 / J3 „ J3 „„^ „ u ^ X 2 (p,q)+p 2 q 2 -X(p,q){p 2 + q 2 +m 2 ) ^ ( (p - qf + : 



J 4 P qE{p)E(q)X(p,q) \{p + q) 2 + mj 

2 f~a -3 , , U ^ X 2 (p, q) + p 2 q 2 + 3X(p,p)q 2 + 3X(q, q)p 2 { (p ~ q) 2 + m 2 , 

~ 9 W 3 pd qg{q)g ] {p) , . , , , r In 

J <ipqE(p)E(q)X(p,q) \{p + q) 2 + m 2 u 

From now on, we use only the MFA and neglect the Fock term. 
The superfluid gap can be determined by using Eqs. ^ and §3§ 

A 3 (p) = + 



E{p)E{q) 16 J \(p + q) 2 +m 2 J pqE(p)E(q) 

S& f*3 ^ , f(P - ?)* + m l \ 2E (P) E (l) ~ m2 



-fj **«*<«> IrMJ fmm (26) 

Finally it must be stressed that in applying Eq. ([3]) for evaluating the quasi- 
particle spectra, the effective mass M* must be considered a functional of the 
distribution function /. 



2.3 Effects of the spatial extension of nucleons 

Standard field theory considers physical particles as structureless, point-like 
objects. This could be a serious shortcoming when composed states in a dense 
medium are described. It is well known that, for instance, the energy and 
charge density of a soliton spreads over a finite range of space [33] . Therefore it 
is legitimate to assign a intrinsical length scale to nucleons immersed in a dense 
environment. 

This was the argument supporting many phenomenological studies of the nuclear 
interaction [551 12S1 1201 OH] ■ From a practical point of view, one can consider N a 
fermions of class a distributed over a finite volume V, then the available space for 
quantization is V = V— ^2 a N a v a , where v a is the spatial extension of this state. 
But the canonical procedure uses the full volume V, this situation is corrected 



8 



by introducing a correction factor y/V'/V in the second quantization of fields. 
This, in turn, modifies fermion bilinears like particle number n a =< ^ a ^a > 
and normal energy density by a factor 9 = V /V = 1 — J^b n v b- The superffuid 
energy density requires a correction O 2 , which is transferred to the formulae of 
the gap function. 

The effective volume inaccessible for other particles due to the presence of a 
spherical object of radius R a is 

V a = U — R a , 

The parameter R a can be understood in a simple minded model as the geomet- 
rical size of a particle. Actually it introduces into the model a characteristic 
scale of the strong interaction, i. e. the spatial spreading of a bounded state 
of quarks and gluons. The value selected for this radius must be compatible 
with similar lengths adopted in hybrid models of the nuclear interaction, see for 
instance [31]. Returning to the schematic picture, the factor a takes account 
of the fact that the inaccessible volume exceeds the actual size of each particle 
and it depends on the spatial arrangement adopted by the collection of objects. 
The minimal volume configuration for identical particles corresponds to a face 
centered cubic arrangement. In such a case is a = 3^/2/%, which is the value 
adopted in the present calculations. 

It is worthwhile to mention that the normalization of the nucleon held with 
an excluded volume coefficient is not equivalent to the introduction of a hard- 
core potential. The normalized nucleon held interacts dynamically with the 
meson fields, both scalar and vector. The in-medium properties of protons and 
neutrons, as well as the meson fields configuration arise simultaneously from 
this interaction. The sigma meson gives rise to the attractive channel of the 
nuclear force, whereas the omega meson is responsible for the repulsive compo- 
nent. Therefore, the proposed normalization affects both attractive and repul- 
sive channels. This fact is evident from Eq. (|9]) and the discussion given above. 
It should be clear that the mean-field value of the sigma meson is strongly af- 
fected by the normalization of the nucleons. Moreover, the relation is highly 
non-linear. Furthermore, the treatment of the mesons is not symmetrical since 
the omega meson is coupled to a conserved charge. In consequence the omega 
meson mean field value is completely determined by the conserved baryonic 
density. The sigma meson mean field value, instead, come forth the hadronic 
dynamics. 

Since introduces an explicit dependence upon the baryonic densities, the 
quasi-particle energy e, see Eq. (J3J, gets an extra term 

2 tal 



2«a y-s V f^V) 



e a (p) = VP 2 + M* 2 + 9^ + ^ M 3 <Z 



3 J H E a {q) ' 

it must be noticed that the additional term depends on density and temperature, 
but not on the momentum. 
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R[im] 


\FEC / ANFEC 
max I max 


9w/g s 


0.6 


0.68 


1.222 


0.62 


0.64 


1.209 


0.64 


0.61 


1.193 


0.66 


0.59 


1.171 


0.68 


0.58 


1.141 


0.70 


0.56 


1.095 


0.71 


0.57 


1.062 


0.72 


0.60 


1.018 



Table 1: The gap function on the Fermi surface and the quotient of coupling 
constants for several values of the characteristic nucleon size R. 

3 Results and Discussion 

The model of nucleons and mesons with finite volume corrections has several 
parameters, we used M = 940 MeV for the mass of the degenerate nucleons, and 
m w = 783 MeV for the omega-meson mass. The sigma-meson mass has been 
fixed at m s = 520 MeV, in agreement with [5]. A discussion about variation of 
m s can be found in [5]. 

The length scale R a is not determined by the model, therefore we consider it as 
a constant value ranging between 0.5 fm and 0.9 fm. 

The coupling constants g s , g w are fixed in order to reproduce the saturation 
properties of symmetric nuclear matter, the binding energy Eb = — 16 MeV 
and a saturation density corresponding to the Fermi momentum pp = 1.42 fm^ 1 
IP- 

To have a look of the performance of QHD models in describing the nuclear 
superfluidity, we have examined two different models of the nucleon-meson in- 
teraction. In first place we study the density dependent coupling model of 
|19j . which translate Dirac-Bruckner calculations with Bonn A potential into 
the covariant field theory. We have used the rational function parametriza- 
tion given there for the couplings, obtaining a maximum value A max (pp ) ~ 15 
MeV at n/riQ — 0.25. As another QHD example we take the non-linear meson- 
nucleon model of [35] . it predicts a lower compressibility and higher effective 
mass in the MFA than the a — ui model jTTJ [TB] does. In our calculations we get 
A max (p F ) = 10.6 MeV at n/n ~ 0.1 . 

As the next step we consider the effects of the finite extension of nucleons, 
taking R as an adjustable parameter. The results are summarized in Table 
1, where the treatments regarding the finite extension of nucleons (FEC) or 
neglecting it (NFEC) have been distinguished. There is a sensible reduction of 
A max , between 30% and 40 %, respect to the previous calculations. A non- 
monotonous dependence on R is obtained, the lowest value A max = 5.5 MeV is 
reached for R = 0.7 fm at pp — 0.8 fm -1 . 

Choosing R = 0.7 fm -1 for subsequent calculations, we compare the momentum 
dependence of the gap function for a fixed density n/rio = 0.25. Results are 
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shown in Fig.l. The square points show the magnitude of the gap at the Fermi 
surface. A comparison of the two cases NFEC and FEC, shows that the latter 
yields the lowest absolute value for q < 1.6 GeV. A reduction of almost 40% is 
registered at the Fermi surface in the FEC case. 
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2 4 6 8 10 
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Figure 1: The gap energy as a function of the momentum, square symbols 
indicate its value at the Fermi surface. 

In Fig. [2] we display A(pp) as a function of the density, here one can see 
that the FEC treatment not only lowers A max , but also reduces more than 20% 
the range of densities where the gap is effective. In this case pp = 1.6 fm -1 is 
the upmost value for non-zero pairing gap. 

To improve understanding of this outcome, we investigate separately contribu- 
tions of scalar and vector character to the gap function. The attractive (v s ) and 
repulsive (v w ) potentials evaluated at the Fermi surface 



v a (q) = - 

v w {q) = - 

where Ef = 
simplified to 



8E F E(q) 



1 



4M 2 - 2{E F - E(q)f 



gj 2E F E(g) - M 2 
8 PF qE F E{q) 



In 



2p F q 
(pf - g) 2 + mj 
(pf + q) 2 + m 2 w 



In 



(Pf - q) 2 + ml 
(pf + q) 2 + m 2 



\Jp 2 F + M 2 , have been defined in order that Eq. ([26|) can be 



A 3 (p; 



dqgs{q) [v s {q) +v w (q)} 



As can be seen in figure [3l both v s and v w appears diminished in the FEC 
case, with a stronger suppression of v s as compared to v w in the high momenta 
regime. Therefore v s + v w asymptotically goes to zero faster in FEC than in 
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Figure 2: The energy gap evaluated on the Fermi surface in terms of the particle 
density. Results including Finite Extension Corrections (FEC) or not (NFEC). 

NEFC, minimizing contributions from higher momenta. 

On the other hand, for q < 1 fm _1 the attractive potential is stronger than the 
repulsive one in FEC than in NFEC approaches. 

There is a surprising numerical similitude between the values described above 
for A max (pF) and some findings in [7], where a sophisticated evaluation of in 
medium effects over the 1 Sq pairing is presented, within the Brueckner theory 
of nuclear matter. The BCS gap function is found to have a maximum value 
A. rnax ~ 5 Mev in symmetric nuclear matter at pf = 0.9fm _1 . The upper limit 
for the existence of the gap is given by pf = 1-5 frn . 

Another interesting comparison can be made with the results found by Mat- 
suzaki et al. |16j . That work, as the present approach, evaluate the possibility 
of giving a unified description of both particle-particle and particle-hole chan- 
nels. In that case the sigma-omega model is used in the mean field approxima- 
tion, and the results are rendered physically acceptable by a direct intervention 
over the integrals in momentum space. This is achieved by the introduction 
of a form factor depending on a single parameter, which is adjusted to obtain 
the best fit to microscopic calculations. It must be pointed out that the cor- 
rection factor used in our procedure can not be rigorously considered as a 
form factor, since for a given Fermi momentum it reduces to a constant value 
and therefore it does not modify the integrals. Despite the procedural differ- 
ences, our results are comparable to those of [16] . For instance, if we consider 
the momentum dependence of the single-particle potential v(pp,q) = v s + v w , 
evaluated at pf such that a maximum of A max (pF) is obtained, it takes values 
-4 MeV < v(pF,q) < 3 MeV in 16 . In our calculations, instead, we found -5 
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Figure 3: The pairing potential and its 
terms of the particle density, for FEC 
approximations . 
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MeV < v(pp,q) < 1 MeV, see Fig. 3. Furthermore, the asymptotic behavior 
is similar in both calculations. More appreciable differences are found for the 
momentum dependence of the gap corresponding to the same pp . At very low 
momentum we have A ~ 6 MeV, then as the transfer momentum q is increased, 
the gap decreases smoothly, passes through zero at q ~ 3.5 fm -1 , reaches a 
minimum value of -1 MeV at q ~ 5.5 fm -1 , and finally tends asymptotically 
to zero from negative values, see Fig. 1. As can be seen in Fig. 3b of [15] . 
the low momentum gap is sensibly lower A ~ 4 MeV, the gap passes trough 
zero at a lower value q ~ 2 fm -1 and reaches the same minimum value but at 
q ~ 3.5 fm . From this observations, we can conclude that in our approach, 
the contributions coming from momenta q < 3 MeV are overestimated as com- 
pared with the treatment of [16] . Consequently, the behavior of A max (pp) has 
a maximum value that exceeds by 2.5 MeV the results shown in Fig. 2a of [15] . 
However, this maximum value is reached for pp ~ 0.8 fm in both cases. The 
spreading is also similar, a drop of about 75% is verified at pp ~ 1.2 fm -1 in 
both calculations. 

As a last application we study the temperature behavior of the gap function 
evaluated at the Fermi momentum in the FEC approach. In fig. [J] we select 
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some definite values of the particle density such that the gap has magnitude 
higher than 0.1 MeV at zero temperature. For the lower densities a steep fall is 
registered around T ~ 2.5 MeV. 
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Figure 4: The gap energy for the 1 Sq phase in terms of the temperature for 
several densities. 



4 Summary 

In this work we have studied the superfluid phase in a nuclear environment, 
within the relativistic field theory of hadrons. The Fermi-Landau liquid frame- 
work, as stated by [35], has been used in order to obtain expressions for the 1 S'o 
superfluid gap energy and the distribution functions for both normal and su- 
perfluid phases. This is equivalent to a Hartree approximation, which provides 
a quasi-particle picture of nucleons dressed by the mesonic interaction, solved 
in a self-consistent way with a BCS scheme for the superfluid phase. 
It is well known that BCS approaches in terms of QHD models, produce exces- 
sively large values for the gap in isospin symmetric nuclear matter. In this work 
we explored the possibility to state a easy to handle formalism, capable to re- 
tain the good properties of the QHD theory but modifying the above mentioned 
failure. With this purpose we took the simplest sigma-omega interaction and 
we introduced short range correlations inspired by the finite volume extension 
of nucleons. In our approach an additional model parameter is required, which 
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is related to the characteristic length of the finite size of nucleons. 
We have found a significative reduction, between 30 — 40%, of the 1 So gap in 
symmetric matter. So, the results obtained are closer to non-relativistic es- 
timates. The modifications proposed modify the high momenta performance 
producing a pairing potential which goes to zero faster than in [3]. Further re- 
finements, such as Fock term, additional mesons, non-linear sigma terms, etc., 
could improve the agreement with expected values. 

Finally, the temperature behavior of the superfluid gap obtained is comparable 
with previous estimates @], major differences correspond to the lowest densities. 
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